Light scattering and phase behavior of Lysozyme-PEG mixtures 
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Measurements of liquid-liquid phase transition temperatures (cloud points) of mixtures of a pro- 
tein (lysozyme) and a polymer, poly(ethylene glycol) (PEG) show that the addition of low molec- 
ular weight PEG stabilizes the mixture whereas high molecular weight PEG was destabilizing. 
We demonstrate that this behavior is inconsistent with an entropic depletion interaction between 
lysozyme and PEG and suggest that an energetic attraction between lysozyme and PEG is re- 
sponsible. In order to independently characterize the lysozyme/PEG interactions, light scattering 
experiments on the same mixtures were performed to measure second and third virial coefficients. 
These measurements indicate that PEG induces repulsion between lysozyme molecules, contrary to 
the depletion prediction. Furthermore, it is shown that third virial terms must be included in the 
mixture's free energy in order to qualitatively capture our cloud point and light scattering data. 
The light scattering results were consistent with the cloud point measurements and indicate that 
attractions exist between lysozyme and PEG. 



The addition of small polymers to a dispersion of large 
colloids can result in precipitation or crystallization of the 
colloids. Even when the sole colloid/polymer interaction 
is steric repulsion, an attraction between a pair of col- 
loids is generated by the exclusion of polymer molecules 
from the region between the colloids. This entropic effect 
is known as depletion attraction 0, Q • The water solu- 
ble polymer poly(ethylcnc glycol) (PEG) has been utilized 
extensively to induce protein crystallization Q . Can the 
mechanism of PEG induced crystallization be explained 
by the purely entropic depletion effect? This work ad- 
dresses the related question of whether or not PEG in- 
duces attraction between lysozyme molecules. 

In addition to a liquid to crystal transition, protein 
solutions exhibit a metastable liquid-liquid phase transi- 
tion when cooled |J, |5], |6J, [7J, |8j . This phase transition 
temperature is termed the cloud point (T c i ou d) since at 
this temperature a transparent protein solution becomes 
turbid as liquid droplets of high protein concentration 
form in a liquid of lower protein concentration. The ef- 
fect of added PEG on T c i ouc i was studied for several of 
the 7-crystallin proteins by Benedek and co-workers who 
found that T cloud of 7S crystallin increases as PEG is 
added 8] in agreement with depletion attraction whereas 
for 7D crystallin their results |9j indicate a departure 
from pure depletion. Galkin and Vekilov 01 studied 
lysozyme/PEG mixtures and found that the effect of 
PEG molecular weight on T c i ou d depended on the so- 
lution's ionic strength so no clear evaluation of the de- 
pletion effect could be made. Because the cloud point 
depends on salt type and concentration 0, 0] , all our 
measurements were performed in the same solution con- 
ditions. 

We model the thermodynamics of a protein/polymer 
mixture by expanding the excess Gibbs free energy (G) of 
a two component solution relative to that of the solvent 
in powers of the densities of the two independent solutes, 
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In Eq.QJ pi = Ni/V [Volume -1 ]; ks is Boltzmann's con- 
stant; T is the absolute temperature; Bij [Volume] are 
the second virial coefficients and [Volume 2 ] are the 
third virial coefficients. Virial coefficients are related to 
integrals of the potential of mean force between molecules 
and depend on T [liij . 

A bidisperse hard sphere mixture was used as the ref- 
erence system for lysozyme/PEG mixtures because hard 
sphere systems provide a natural scale for virial coeffi- 
cients. Hard spheres cannot interpenetrate, but have no 
other interactions. The virial coefficients in Eq.Q for a 
bidisperse hard sphere mixture are |l4) : 
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where ri are the hard sphere radii. 

The ratio of lysozyme to PEG measured hydrody- 
namic radii, rjj, in our experiments varied in the range 
3 > ^h S / 7 'h EG — 0-8- Therefore we need to account 
for the fact that PEG molecules, which are nearly as 
large or larger than lysozyme molecules, are not spheri- 
cal and wrap partially around proteins thereby reducing 
the polymer/protein excluded volume. This is done by 
defining an ideal effective polymer radius (r e s) which is 
smaller than the polymer's radius of gyration r g . Eiscn- 
riegler et al. |15| found a closed formula for computing 
r e ff from the protein radius and r g (see Fig. 6 in ref. 15]). 
Our procedure was to equate the protein's radius with its 
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measured rjj whereas the process to determine the PEG's (R2) can be written as: 



r g is described in the following paragraph. 

The PEG radii of gyration for the molecular weights 
employed here were too small to measure with static 
light scattering because r g is much less than the wave- 
length of light. In order to determine r g , we therefore 
first measured r# and then used the relation r g /rH = 
1.48M 2 - 012 [HI where M 2 [g mol" 1 ] is the PEG molecu- 
lar weight. This should be compared to the theoretical 
value, r g /rn = 1.56, for a polymer in a good solvent 
[l7|. For lysozyme we found: r^ s = 2.2 nm which falls 
within the range of previously reported values . We 
measured r^ G =0.75 and 2.7 nm for PEGlk and 8k re- 
spectively where the manufacturer's stated value for M2 
is designated with the nomenclature PEGnk, meaning 
PEG nxlO 3 [g mol -1 ]. The procedure outlined above 
yields r cS =1.2 and 3.8 nm for PEGlk and 8k. 

We measured the cloud point temperature of lysozyme 
solutions as an approximation of the spinodal decompo- 
sition temperature since the cloud point closely tracks 
the spinodal • Starting from the free energy, we calcu- 
late how the spinodal temperature changes with added 
polymer concentration at fixed protein concentration. In 
all that follows the subscript 1 refers to lysozyme and 
the subscript 2 refers to PEG. A two component mixture 
at fixed concentration undergoes spinodal decomposition 
when the temperature reaches T sp defined by [Tjfl : 
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Imposing the constraints that the solution remains on 
the spinodal curve when polymer is added and that the 
protein concentration is constant one finds dT sp /dp2 = 
-(df/dp 2 )/(df/dT). From Eq.JTJ}, one then obtains to 
first order in p\\ 



lim^ 

P2^0 dp2 



— 3Cn2 + 2(£?i2 + 3piCn2) 



dB n /8T+3 Pl 



dT 



(4) 



In all that follows we have set dCm/dT = since ex- 
perimentally Cm = at all temperatures. 

We define the numerator of Eq.(gJ to be 7. The PEG 
molecular weight dependence enters Eq. only through 
the virial coefficients. From Eqs.(|21IBlEJl, it can be shown 
that for a hard sphere mixture 7 > for all sphere sizes 
7'i and T2 and concentrations p\. In the limit of small 
polymers (q = r 2 /ri < 1) we find that 2-f / (B^ / 4) 2 = 
12g 3 in agreement with scaled particle depletion models 

In order to explore the dependence of T c i ou d on PEG 
molecular weight further, virial coefficients were deter- 
mined from light scattering experiments on the same 
lysozyme/PEG mixtures and then compared with pre- 
dictions of the depletion theory. Kirkwood & Goldberg 
[20| showed that the excess light scattering of a two so- 
lute system (.R1+2) over that of a single solute system 
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Here K = 2(wn ni) 2 /NaX 4 where n is the solvent re- 
fractive index, m = 4^ is the refractive index increment 
of solute i, Na is Avogadro's number, A is the wavelength 
of the incident radiation in vacuum, c[g mL -1 ] is the so- 
lute weight concentration and R is the Rayleigh ratio. 

a and (3 depend on the added polymer properties (M2 , 
n 2 ), concentration C2, and the protein/polymer interac- 
tion El: 
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B12 is obtained from measurements of a as a function of 
polymer concentration C2 shown in FigEI 
The coefficient j3 is given by 
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C112 is obtained from measurements of (3 as a function 
of C2 , as shown in Fig. |5J since all the other quantities in 
Eq.Q are determined independently. 

It is possible to view the two component poly- 
mer/protein solution as an effective one component pro- 
tein solution. The effective protein/protein second virial 
coefficient, -Bjf, can be obtained from Eq.Q by imagin- 
ing the addition of invisible polymers (ri2 = 0) to the 
protein solution [U H3 which yields: Bff = B n + 
c-2 [(3Cn2 — 2£> 2 2 M2)] . Experimentally we do not index 
match the polymer (ri2 7^ 0). Instead we measure the 
virial coefficients and calculate B\i . 

Hen egg white lysozyme was purchased from Seika- 
gaku America. Mi based on sequence is 14,400 g mol -1 . 
Poly(ethylene glycol) (PEG) was purchased from Sigma 
and Fluka. 

The protein and PEG were dissolved in a 0.2M sodium 
phosphate buffer with NaCl 0.5M at pH 6.2 where M = 
[mol L -1 ]. All solutions were centrifuged at ~ 12000 x g 
for 1 hour and then passed through 0.2/im filters directly 
into precleaned scattering cuvettes. The lysozyme con- 
centrations were measured by UV absorption using an 
extinction coefficient £280nm — 2.64 mL mg" 1 cm" 1 . 

Refractive index increments of lysozyme and PEG were 
measured using a Brookhaven Instruments differential 
refractometer at A = 620 nm. For lysozyme dn/dc — 
1.85 x 10~ 4 mL mg- 1 and for PEG dn/dc = 1.34 x 10~ 4 
mL mg _1 independent of M2 and T. 

Cloud point temperatures were determined by opti- 
cal microscopy. Rectangular glass capillaries (0.1 mm 
path length, VitroCom) were filled with solution, flame 
sealed and then placed in a custom built temperature 
controlled microscope stage. The temperature at which 
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FIG. 1: Cloud point temperatures (T c i ouc j) of lysozyme/PEG 
mixtures as a function of PEG concentration for different 
PEG molecular weights ranging from 400 to 35k g/mol at a 
fixed lysozyme concentration of 48 mg/mL in the phosphate 
buffer. The lines are guides to the eye. 



homogeneous nucleation of dense droplets occurred was 
called T c i ou( j. The temperature of each solution was cy- 
cled up and down at approximately 1°C min -1 through 
Tcioud severa l times for each measurement. Little or no 
difference (< 0.5°C) was observed when comparing T c loud 
obtained by cooling and heating. 

The static and dynamic light scattering experiments 
(SLS & DLS) were performed using an ALV goniometer 
and correlator system in the vu polarization mode. Abso- 
lute Raylcigh ratios of aqueous solutions were determined 
by using pure toluene as a standard whose Rayleigh ra- 
tio is known |2^] . The lysozyme and PEG hydrodynamic 
radii th were obtained from DLS measurements on dilute 
solutions 

Our T c i ou d measurements of lysozyme/PEG mixtures 
at constant ionic strength reveal a systematic trend: 
whereas T c ioud increases upon the addition of high 
molecular weight PEG, it decreases for low molecu- 
lar weight PEG as shown in Figure ^ We extracted 
lim P2 ^o dT c \ ou d/dp2 from this data. Measurements of 
Tcioud made at higher lysozyme concentrations near 
lysozyme's critical point showed that \im P2 ^odT/dp2 is 
independent of lysozyme concentration (data not shown) . 
The hard sphere mixture model predicts that the cloud 
point temperature dependence on PEG concentration 
cannot change sign with PEG molecular weight, i.e. 
7 > 0. However, our experiments reveal qualitatively 
different behavior. For dilute PEG concentrations, low 
molecular weight PEGs depress the cloud point, i.e. sta- 
bilize the solution, which is opposite to the depletion pre- 
diction, while higher molecular weight PEGs raise the 
cloud point, i.e. destabilize the solution. 

Table [I] shows the PEG molecular weights obtained 
from SLS data. PEGlk and 8k second virial coeffi- 
cients as functions of temperature were found to be 
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FIG. 2: The variation of a (Eq.©) and (3 (Eq.Q) on PEG 
concentration (C2) at T=30°C is shown for lysozyme/PEGlk 
and PEG8k mixtures dissolved in the phosphate buffer. The 
solid lines in panel a) indicate the fits to a used to obtain B12 
by Eq.lJfjJ. The solid lines in panel b) indicate the fits to f3 
used to obtain C112 by Eq.Q. The dashed lines in panel b) 
show the dependence of /3 on PEG concentration assuming 
Ciia = 0. 



Sf 2 lk (T)[mL mol g~ 2 ] = 0.01 - (8.8 x 10~ 5 x T[°C]) and 
Sf 2 8k (T)[mL mol g~ 2 ] = 5.6xl0- 3 -(9.5xl0- 5 xT[ o C]), 
which agree well with previous results [3. Tabledshows 
that measured PEG second virial coefficients are the 
same order as those for equivalent hard spheres, imply- 
ing that PEG interactions are repulsive. For lysozyme: 
Bu(T = 20°C) = -3.95 ± 0.3 x 10~ 4 mL mol g" 1 , 
dBn/dT = 1 x 10" 5 mL mol g~ 2 °C" 1 and M x = 
13800 ± 500 gmor 1 . 

Eq. JfjJ) demonstrates that important information about 
lysozyme/PEG interactions is contained in a[c2), namely 
B12. The measured values of B12 shown in Table Q]are 
consistent with previous measurements [25L . If inter- 
ference between the scattering from lysozyme and PEG 
was negligible then da/dc2 = 0, as seen from Eq.JSJ) 
with (M2fi2) / (Mini) = 0. However, Figure |3 panel a) 
shows that da/dc2 > consistent with Prausnitz and co- 
workers |25| . This demonstrates that interference from 
lysozyme and PEG scattering may not be ignored as done 
by Kulkarni et al. |27l | and that a may not be treated as 
a constant. Figure |2] panel b) displays the variation of 
(3 with PEG concentration. The solid lines through the 
data points are the linear fits used to obtain the values 
of C112 by Eq.© shown in Table H We found B X2 and 
C112 to be temperature independent. 

Much interest has been generated by the conjecture 
that the second virial coefficient may be sufficient to pre- 
dict protein solution phase behavior [28l Ej| . The mea- 
sured dependencies of T c i olK j and (3 on PEG concentration 
show that third virial coefficients must be included in the 
free energy of lysozyme/PEG mixtures, i.e. C112 7^ 0. 
If C112 = then Eq.0 predicts (3{c2) shown as the 
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dashed lines in panel b) of Fig0 which disagree with 
our data. Additionally, if Cm = then Eq.QJ predicts 
that dT sp /dp2 cannot change sign as a function of PEG 
molecular weight, which is inconsistent with Fig^ 

Table HI compares our measured values of the mixed 
virial coefficients, B i2 and Cn 2 , to Bf 2 s and C^l, cal- 
culated from Eqs.(0El)- The lysozyme equivalent hard 
sphere radius was taken to be rjj and the PEG equiv- 
alent hard sphere radii were taken to be r e g- The 
lysozyme/PEG mixed virial coefficients are less than 
those of equivalent hard sphere mixtures. Therefore, at- 
tractive interactions must exist between the protein and 
polymer or the repulsion must be less than that between 
hard spheres. 

The measured virial coefficients allow us to calculate 
the variation of the effective lysozyme second virial coef- 
ficient with PEG concentration, dBlf/dc 2 . Depletion 
theory predicts that PEG induces attraction between 
lysozyme molecules in which case dBif/dc 2 < as seen 
from Eqs.0EJ. Conversely, if &Bff/dc 2 > then PEG 
induces repulsion between lysozyme molecules. For PEG 
8k we find dB\f/dc 2 = 0.33 ± 0.2 x 10~ 3 mL 2 mol g~ 3 
and for PEG lk we find that dBf^/dc 2 = 1.6±0.2 x 10~ 3 
mL 2 mol g~ 3 . We conclude that adding PEG weakens 
the attraction between lysozyme molecules in contradic- 
tion to depletion theory. 

We find experimentally lim C2 _ 9T c i ou d/<9c2 = -0.15, 
0.06 ±0.02 °C mg- 1 mL- 1 for PEG lk and 8k re- 
spectively. Using Eq.Q with the measured virial 
coefficients obtained independently from SLS yields 
lim O2 ^ dT sp /dc 2 = -0.03±0.05, 0.005±0.008 °C mg- 1 t 4 
mL -1 for PEG lk and 8k respectively. These results 
show that the measured virial coefficients can correctly 
reproduce the sign of lirm, 2 _»o <9T c i ou d /dp 2 ■ As discussed 
previously, 7 > for a hard sphere mixture whereas Fig^ 
shows that experimentally the sign of 7 varies with PEG 
molecular weight. Therefore a depletion model cannot 
account for the observed variation of the cloud point on 
polymer concentration, but the measured virial coeffi- 
cients in conjunction with Eq.Q do so qualitatively. 

This paper demonstrates that the depletion theory 
does not describe PEG/lysozyme mixtures. Firstly, the 
observation that adding low molecular weight PEG de- 
presses T c i ou d whereas high molecular weight PEG raises 
Tcioud cannot be accounted for by a pure depletion model. 
Secondly, depletion theory predicts that adding PEG in- 
duces an attraction between lysozyme molecules whereas 
light scattering revealed the opposite: PEG induces re- 
pulsion between lysozyme molecules. Lysozyme/PEG 
interactions were characterized by virial coefficients ob- 
tained from light scattering experiments. It is demon- 
strated that to even qualitatively explain our data, the 
free energy must include third virial terms. The mea- 
sured mixed virial coefficients are smaller than those pre- 
dicted for an equivalent hard sphere mixture and are 
consistent with attractions between lysozyme and PEG. 



Models of hydrogen bonding of water molecules to the 
PEG backbone qualitatively explain PEG's phase be- 
havior alone in water j^. Therefore we speculate that 
PEG can similarly form hydrogen bonds with residues on 
the surface of lysozyme molecules thereby creating the 
attraction between PEG and lysozyme. The measured 
virial coefficients, combined with thermodynamic theory 
predict the observed behavior of the cloud point demon- 
strating the consistency of these two independent sets of 
experiments. Therefore, an accurate model of the phase 
behavior and light scattering of lysozyme/PEG mixtures 
must account for both the entropic depletion effect and 
an energetic attraction between protein and polymer. 
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Mi 
[10 3 g/mol] 


N 


B22 

[10 _3 mL mol/g 2 ] 


B22/BW 


Bi 2 

[10" 4 mL mol/g 2 ] 


b 12 /bW 


C112 

[10" 4 mL 2 mol/g 3 ] 


C112/CU2 


PEGlk 


0.97 ± 0.1 


22 


9.2± 1.0 


0.5 


12.9±1.2 


0.34 


16.1±1.3 


0.23 


PEG8k 


10.4±0.5 


236 


3.35±0.12 


0.4 


1.7±0.3 


0.089 


3.13±0.2 


0.051 



TABLE I: PEG molecular weights (M2) and second virial coefficients (B22) taken from light scattering data at T= 30°C, the 
degree of PEG polymerization (N = A/2/44) and mixed virial coefficients (Bi2,Cii2) for lysozyme/PEG solutions obtained 
from fits to a and f3, Eqs, l6l7ll . are shown. The equivalent hard sphere values were obtained by using r g of PEG to determine 
B22 an d r e s of PEG along with rjj of lysozyme to determine B 12 and Cn% by Eqs, l2l3ll . The error bars came from the linear 
least squares fits. 
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